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Abstract. We study the decoherence of atomic interferometers due to the scattering of stochastic gravita- 
tional waves. We evaluate the 'direct' gravitational effect registered by the phase of the matter waves as 
well as the 'indirect' effect registered by the light waves used as beam-splitters and mirrors for the matter 
waves. Considering as an example the space project HYPER, we show that both effects are negligible for 
the presently studied interferometers. 

PACS. 03.65.Yz Decoherence; open systems; quantum statistical methods - 03.75.-b Matter waves - 
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1 Introduction 

The idea that spacetime fluctuations could play a uni- 
versal role in the transition from quantum to classical 
physics has been proposed by a number of authors. Al- 
ready present in the Feynman lectures on gravitation 1 , 
it was more thoroughly developed and popularized for 
instance in |§,§J^,^|. An important argument in favor of 
such an idea is that the Planck mass, i.e. the mass scale 



mp = J!% built on the Planck constant H, the velocity 
of light c and the Newton constant G, has a value mp ~ 
22 [ig lying on the borderland between microscopic and 
macroscopic masses. In other words, microscopic masses 
could be characterized as masses m < mp for which the 
associated Compton length Iq — — is larger than the 

Planck length lp — J ^ ~ 10~ 35 m whereas macroscopic 

masses m > mp would correspond to a Compton length 
£c smaller than the Planck length £p. 

Clearly, this dimensional argument is not by itself suf- 
ficient to reach definitive conclusions. It can be hoped 
that the existence of fundamental spacetime fluctuations, 
with a length scale determined by tp, may be revealed by 
long-term diffusion effects in the same manner as micro- 
scopic molecular motion is revealed by Brownian motion. 
It has for instance been proposed that intrinsic spacetime 
fluctuations could be observed through a decoherence ef- 
fect which mieht be visible with matter-wave interferom- 
eters ]^,|^,p|, pi3|]Tl| ] . The effect has not been seen in exist- 
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ing matter- wave interferometers [ |12[jl3| , |14| . More sensitive 
instruments are being developed, like the atomic interfer- 
ometer HYPER designed to measure the Lense-Thirring 
effect in a space-borne experiment, and it is important 
to estimate the ultimate decoherence due to fundamental 
spacetime fluctuations for such an instrument JD|. 

The aim of the present paper is to give quantitative 
answers to this question by considering the decoherence 
mechanism associated with the scattering of gravitational 
waves present in our celestial environment. These gravi- 
tational waves are the intrinsic fluctuations of spacetime 
predicted by general relativity Jl~6|] . The latter theory can 
be used as an accurate effective theory of gravitation for 
all frequencies ever explored in experiments [jl~7| , p~8]| . The 
gravitational waves are the intrinsic field fluctuations pre- 
dicted by the linearized form of the theory p9| , ^0| , pT| . 
This linearized form is widely used for studying propaga- 
tion of gravitational waves and their interaction with the 
presently developed interferometric detectors |p2| , p3|j2^| . 

In the present paper, we will study the decoherence of 
an atomic interferometer due to its interaction with the 
stochastic background of gravitational waves emitte d by 
unresolved sources in our galaxy or its vicinity |^5| , p6| . 
These waves dephase differently the matter-waves on the 
two interfering paths as well as the light-waves used to 
build up beam-splitters and mirrors in atomic interfer- 
ometers [^7j . When averaged over the integration time of 
the measurement, the differential dephasing results in a 
loss of contrast of the interference fringes. We will use the 
remaining fringe contrast to characterize the decoherence 
and write it in terms of the geometry of the interferometer 
and of the statistical function describing the gravitational 
environnement . 
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In this approach, decoherence will be understood as 
resulting from a phase dispersion due to the unobserved 
degrees of freedom of the gravitational environnemcnt. To 
be more precise, gravitational waves with frequencies in 
the detection window of the interferometer have to be in- 
terpreted as signals while frequencies outside the detec- 
tion window are ignored. In a typical situation, the de- 
tection window corresponds to frequencies smaller than 
the inverse of an averaging time. The integration over fre- 
quencies outside the detection window can thus be iden- 
tified with the trace over the environnemental degrees of 
freedom usually considered in theoretical studies of de- 
coherence [|§(29||^,|I|,||||,||,||. The phase dispersion 
approach used in the present paper is known to be equiv- 
alent to the other approaches to decoherence and it is 
obviously well adapted to the description of interferome- 
ters where the phase is the natural variable. 

We will show in this paper that the scattering of grav- 
itational background does not lead to an appreciable de- 
coherence effect for the atomic interferometers presently 
studied, HYPER being chosen as the typical example. In- 
cidentally, this means that atomic interferometers will not 
have their interference fringes destroyed by this decoher- 
ence mechanism. This answer has to be contrasted with 
recently published results which prove that the scatter- 
ing of stochastic gravitational waves present in our galac- 
tic environment is the dominant, and extremely efficient, 
decoherence mechanism for macroscopic motions, say the 
planetary motion of the Moon around the Earth |3^j38[ . 
This contrast is easily explained by the already evoked 
dimensional argument : gravitational decoherence effects 
are likely to be more efficient for macroscopic masses than 
for microscopic objects. In particular, the mass of the 
Moon is larger than Planck mass by orders of magnitude 
whereas the microscopic entities used as spacetime probes 
in atomic interferometers have their mass much smaller 
than Planck mass. However, as already stated, this simple 
scaling argument is not by itself sufficient to answer quan- 
titative questions about the decoherence rates. In the fol- 
lowing, we will give precise estimations of the decoherence 
effect which depend not only on the mass of the atoms, 
but also on their velocity, on the geometry of the inter- 
ferometer and on the noise spectrum characterizing the 
gravitational background in the relevant frequency range. 



In the present paper, we will adopt the common strat- 
egy of studies of gravitational waves : admitting that this 
point can be dealt with, we then chose a specific gauge, 
namely the transverse traceless (TT) gauge with metric 
perturbations differing from zero only for purely spatial 
components hij = I , 2, 3 stand for the spatial indices 
whereas represents the temporal index). Then gravita- 
tional waves are conveniently described through a mode 
decomposition |3§] : 



hij( X ) 



d 4 fc 



h-oo = hoi = 



hij[k]e 



(1) 



Any Fourier component is a sum over the two circular 
polarizations : 

Gravitational waves correspond to wavevectors k lying on 
the light cone and they are transverse with respect to this 
wavevector : 







c 



(3) 



The gravitational polarization tensors are obtained as prod- 
ucts of the polarization vectors well-known from elec- 
tromagnetic theory. When necessary, we will chose the 
following representation for the unit vector n along the 
propagation direction of the gravitational wave and the 
corresponding polarization vectors : 



n 



ck 

U! 



sin 9 cos (f 
sin 9 sin ip 
cos 9 



uj>0 



e' K = 



cos 9 cos ip + sin ip 
cos 9 sin ip — ij cos ip 
sin 9 



(4) 



Spatial vectors are written as bold letters. Note that real- 
ity conditions for the perturbation metric hij (x) are read : 



2 Gravitational backgrounds 

A first step is to characterize the fundamental fluctuations 
of spacetime and their effect on the motion of matter. 
Although a complete quantum theory of gravity is not 
available, it is possible to describe spacetime fluctuations 
in our environment. At the frequencies of experimental 
interest, which are much smaller than Planck frequency, 
they are identified |57]j38| as the stochastic backgrounds 
of gravitational waves currently studied in relation with 
the development of gravitational wave detectors [^] . 

The effect of gravitational perturbations may in princi- 
ple be described in a manifestly gauge-invariant manner. 



l-k]Y 



= W\ 



<[k] 
-k] 



(5) 



We use the natural caracterization of the stochastic 
background in terms of the spectral density of strain fluc- 
tuations Chh[k] of the metric. For simplicity, we consider 
the case of gaussian, stationary, unpolarized and isotropic 
backgrounds : 

h?[k]W'[k']} = (2n)H~ n ' 5 4 (k + k')C hh [k] (6) 

The general case could be dealt with by considering arbi- 
trary correlations between amplitudes h + [k] and h~ [k]. 
This would allow one to take into account polarized or 
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anisotropic backgrounds, as is necessary for a thorough 
analysis of the galactic background, as well as non station- 
nary fields appearing in some cosmological models f2Fj,|JJ . 

Gravitational backgrounds are usually written in terms 
of one metric component (say hi2, but the result would 
be the same for other components due to the isotropy 
assumption) at a fixed spatial position (say x = 0, but 
the result would be the same for other positions due to 
the stationarity assumption) as a function of time t. They 
are thus described by the spectral density 5ViM of strain 
fluctuations considered in most papers on gravitational 
waves detectors [E3l : 



(Mf)Mo)) 



(7) 



This noise spectrum, written in the TT gauge, is not gauge 
invariant. This is not a problem since only gauge invariant 
quantities will be computed in the following. It follows 
from equations (0J^) that Sh [w] is obtained by integrating 
Chh[k] over the momenta k which correspond to a given 
frequency oj : 



S h [u] = J 



d 2 n 



(eJ[k}eJ[k}eJ[-k}e y 2 [-k}C hh [k]) i 



<»n = / m 



d n = dcos 



(8) 



We have denoted by < ■ • ■ > n the averaging over spatial 
directions of wavevectors at a given frequency. Using the 
fact that gravitational wavevectors lie on the light cone 
and the simplifying assumptions already described, we in- 
verse the preceding relation to obtain : 



C hh [k] = l^cH{k 2 )^ 



(9) 



In this paper, we will consider the binary confusion 
background describing gravitational waves emitted by un- 
resolved binary systems in our galaxy or its vicinity. This 
background is represented for example on figure (I) of 
reference It relies on the laws of physics and astro- 
physics as they are known in our local celestial environ- 
ment. In particular, it depends on the statistical reparti- 
tion of binary systems in the sky. In the frequency range 
of interest, which is discussed in more detail below, the 
binary confusion background dominates other sources of 
stochastic gravitational waves, in particular those associ- 
ated with cosmological contributions. Note that this fre- 
quency range, between /iHz and mHz, corresponds to fre- 
quencies much smaller than the detection window of VIRGO 
and other ground-based optical interferometers but roughly 
of the same order as the frequencies in the detection win- 
dow of the space-borne project LISA [JJ 



The properties of the gravitational bath may also be 
characterized by an effective number n gw of gravitons per 
mode (the precise relationships between n gw [u>] and Sh [u>] 
are given in the last section). It is worth noting imme- 
diately that this number n gw is extremely large, so that 
the gravitational environment corresponds to the limit of 
high-temperature classical fluctuations. As a consequence, 
the vacuum fluctuations of the gravitational field Eoj , 
which have been shown to lead to ultimate fluctuations of 
geodesic distances of the order of Planck length |4l| , f42| , [f3| , 
are ignored in the present paper. We also note that, as the 
number of unresolved binary systems contributing to the 
binary confusion background is large and as these sources 
are independent of each other, it appears quite safe to 
consider that the gravitational background obeys gaussian 
statistics (this might not be true for cosmological contri- 
butions). 

In order to discuss our main results in the end of the 
present paper, we will use the fact that the binary con- 
fusion background corresponds to a nearly flat function 
Sh[uj] in the frequency range 1 — lOO^Hz. This means that 
the gravitational noise spectrum is quasi-thermal at such 
frequencies and entails that the decoherence mechanism 
can be interpreted as a Brownian-like diffusion process. 
We will estimate quantitatively the quantity playing the 
role of the diffusion coefficient and show that it does not 
lead to an appreciable decoherence effect for atomic inter- 
ferometers like HYPER. 



3 Gravitational decoherence 

Now we want to quantify the effect of the previously dis- 
cussed gravitational waves on the coherence properties 
of an atomic interferometer. In the present section, we 
present a detailed discussion of this effect which could ba- 
sically be schematized as follows : coherence of the inter- 
ference fringes is preserved if and only if the differential 
phase perturbation between the two arms is controlled to 
a level much better than 2ir. 

When propagating in spacetime, the atomic probe field 
registers curvature fluctuations. The main effect of the 
perturbation is described by the eikonal approximation, 
valid for wavevectors k of the gravitational wave much 
smaller than wavevectors K of the probe field. This effect 
is characterized as a dephasing <I> of the probe field, eval- 
uated at a fixed spatial position (say x = 0) as a function 
of time t. At the lowest order, the dephasing is linear in 
the metric and can be decomposed over the gravitational 
wave modes : 



*(*) 



d 4 fc 

WT 4 



(10) 



7=± 



1 Informations 
to rs may be 



littp:/ /www. virgo.infn.it 



hannover.de^ 



on the ground-based GW deteo 
found on the Web sites; VIRGO 

G EO : tittp://www.gco600.un i- 
LIGO : [rttp://www,ligo.caltcch,cdu 



TAMA : |http'://tamago. mtk, nao, ac.jp/tama, htmT ; ACIGA 



|http://www.anu.cdu.au/Physics/ACIGA^ 



Informations on the space-borne inte rferometer LISA may be 



found on the Web sites at NASA : 


http://lisa.jpl.nasa.gov/ 


and ESA : 


trfctp: / /sci.esa.int /home/lisa'-- 
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At the moment, 4>J.(t) are time-dependent coefficients which 
depend on the geometry of the interferometer. These co- 
efficients will be explicitly written in the forthcoming sec- 
tions but we already know that they satisfy the reality 
conditions : 

<t>i k = m (ii) 

For a stationary, isotropic and unpolarized background, 
the correlation function for the dephasing is deduced : 

(*w*(*o> = / (0 E m4>i k (t')c hh [k]e-^-^ 

(12) 

The gravitationally induced dephasing is not always a sta- 
tionary noise. Here, we will focus our attention on situa- 
tions where it is stationary or quasi-stationary and where 
the correlation function is simply represented by a noise 
spectrum : 

(#(*)*(*')) = / ^#Me- ,u( ^ f,) 

J 27T 

S$[w] = S h [u)]A[w] , ^H = |E(l^l 2 )n (I 3 ) 

7=± 

In fact, we have supposed that the amplitudes 4>1(t) are 
only slowly dependent functions of time and, furthermore, 
we have used the reality condition . The equation thus 
obtained means that the stochastic dephasing has a noise 
spectrum equal to the product of two factors, the gravi- 
tational noise spectrum Sh [oj] and the apparatus response 
function A[uj] which depends on the geometry of the ap- 
paratus. 

As mentioned earlier, the analysis of decoherence should 
in principle take into account the detection strategy used 
in the interferometric measurement. In order to fix ideas, 
we will consider a simple strategy where the output of 
the interferometer, supposed to be linear in the variation 
of the dephasing is averaged over an averaging time 
T av . This means that the signal, for instance the Lense- 
Thirring effect in HYPER experiment |l5| , is contained in 
the averaged dephasing <P which is read in the frequency 
domain as : 

= — ^r- $M , r=— (14) 

I — ILO T av 

In other words, the signal window is defined by a low-pass 
filter with a bandwidth r. Now, the frequencies outside 
the signal window constitute an uncontrolled noise which 
may degrade the coherence of the interferometer if it is 
large enough. 

In order to estimate this potential decoherence effect, 
we define the uncontrolled dephasing which is also the de- 
phasing after a high-pass filter with the same bandwidth : 

8$[u) = <P[lu] - = ~ tW (15) 



We then consider the visibility V of the fringes which is the 
mean value of the exponential of the uncontrolled phase 
noise : 

V= (exp (16) 

As we have supposed the gravitational background to obey 
the gaussian statistics, the various dephasings are also 
gaussian stochastic variables. It follows that the visibility 
V of the fringes may be expressed in terms of the variance 
of the uncontrolled noise : 

V = cxp(-^) , A$ 2 = (5$(t)5$(t)) (17) 

Finally this variance is given by the following integral over 
frequency : 

A * 2 = J ^ Sh[uj]A[io] u^Tn (18) 

The last two equations are a fondamental result of this 
article. They quantitatively characterize the gravitational 
decoherence effect through the reduction of the fringe vis- 
ibility V. This visibility is the exponential of the phase 
noise variance. The latter is an integral over the whole 
frequency spectrum of a product of factors. Besides the 
gravitational noise spectrum Sh [u>] and the apparatus re- 
sponse function A[lo] which have already been discussed, 
there is a further factor, a high-pass filter with a cutoff 
r, which defines uncontrolled noise as corresponding to 
frequencies outside the detection window. For simplicity, 
we have considered here a Lorentzian expression associ- 
ated with a simple averaging strategy. Note that the filter 
may help us to regularize potential infrared divergences 
in the forthcoming calculations. We shall see later on that 
this is not necessary for the situations considered in this 
paper and that the results are essentially independent of 
the cutoff. Note also that it would be easy to replace the 
Lorentzian filter by more sophisticated expressions corre- 
sponding to different signal detection strategies. 

As it is usual in decoherence theory, the remaining co- 
herence, here the visibility V, is the exponential of a noise 
variance, the variance A<P 2 of the uncontrolled noise. From 
the point of view of the interferometrist, this expression 
has a quite clear expression : the gravitational perturba- 
tion dephases differently the waves in the two arms, and 
the fringe contrast is appreciably degraded if and only 
if the resulting variance A& 2 is of the order or greater 
than unity. At this point, it is worth emphasizing that 
the same result would have been obtained through more 
formal approaches of decoherence with the uncontrolled 
noise frequencies thus interpreted as the degrees of free- 
dom of the environment (see for example pq] for a more 
detailed discussion of this point). 

Expression ( |l8|) is valid for a gaussian, stationary, isotropic 
and unpolarized background. The numerical factor | is 
a factor arising from angular averaging for an isotropic 
and unpolarized background. As already mentioned, gen- 
eralization to anisotropic, polarized and non stationary 
background is possible. Even non gaussian noise can in 
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principle be dealt with by developping the fringe visibil- 
ity in cumulants of the stochastic noise. However, these 
refinements would not change the main conclusion of the 
paper, namely that the gravitational decoherence is not 
efficient for presently studied atomic interferometers such 
as HYPER. The next sections are devoted to explicit cal- 
culations of the apparatus function A[u>] for such an in- 
terferometer. 



4 Atomic interferometer with a 
Mach-Zehnder geometry 

In order to smoothly introduce the more technical parts 
of our evaluations, we first consider in this section an 
hypothetical atomic interferometer with the same Mach- 
Zehnder geometry as HYPER but where optical elements, 
beam splitters and mirrors, would be built up from mas- 
sive and motionless material pieces. A more realistic de- 
scription of HYPER with optical elements built up on 
stimulated Raman processes is presented in the next sec- 
tion. 

The Mach-Zehnder geometry is represented on figure 
(|l|). The atomic matter- waves are supposed to follow their 
classical trajectories. Equivalently, they obey the propa- 
gation equation of a scalar field with a wavevector and 
the dephasing <P at is the first-order effect of the gravita- 
tional perturbation of the associated lagrangian. Then, the 
dephasing evaluated on a closed loop is gauge- invariant, 
as soon as conservation of the energy-momentum tensor 
is properly taken into account [Q. This raises specially 
delicate problems for the description of beam splitters and 
mirrors in the interferometer. 

In the present section, we solve these problems by us- 
ing the method commonly adopted for the description of 
optical interferometers [pMq] : the treatment is largely 
simplified when optical elements such as beam splitters 
and mirrors are described as heavy objects initially at 
rest and gravitational waves described in the TT gauge. 
In this gauge indeed, gravitational waves have no effect 
on massive objects at rest. And objects with a large mass 
may be considered as staying at rest when initially at rest, 
since the momentum transfered by the field upon scatter- 
ing does correspond to a negligible velocity change. This 
implies that the effect of gravitational waves on optical 
elements can be ignored throughout the calculation. This 
treatment is certainly an approximated one. In particular, 
the assumptions just discussed entail that the interferom- 
eter is at rest. This is obviously not the case for HYPER 
which is orbiting around the Earth, this motion being im- 
portant for the analysis of the signal associated with the 
looked for Lense-Thirring effect jl5| . However this motion 
may be disregarded in the discussion of decoherence in 
the present paper. In the present section, we build up our 
calculations on this approximation. 

Now we write the dephasing <£ at (t) obtained from the 
geodesic deviation equation for the probe field : 



Ko 
2 



t + T 



hijit'^it'Wit^cdt' (19) 



This equation is written in the eikonal approximation where 
the wavevector fc M of the gravitational waves is much smaller 
than the wavevector of the probe field f47j . The inte- 
gral ( jl9| ) is taken along the unperturbed geodesic path 
x(t') of the probe with the coordinate time t' used as an 
affine parameter; r denotes the time of propagation along 
the path, and u is the reduced wave vector : 



u 



K 



(20) 



We then deduce the coefficients 4>1(t) which characterize 
the sensitivity of the dephasing to the gravitational wave 



m = 



Kg 
2V2 



t+T 



(e 7 [/c]*.u(0) 2 e-**"*^*' W (21) 



These coefficients are time-independant, due to the sta- 
tionarity assumption and to the fact that we have sup- 
posed the interferometer to be at rest. 




Fig. 1. Schematic description of an atomic interferometer with 
a Mach-Zehnder geometry : thanks to the presence of beam 
splitters at A and D and mirrors at B and C, atomic matter- 
waves are coherently recombined after having propagated along 
the two geodesic paths ABD and ACD; the output intensity 
is then a sinusoidal function of the differential dephasing # 
between the two paths. For simplicity, we consider that the 
interferometer has the symmetry of the rhomb. The angle 2a 
is exaggerated on the figure. 



As a consequence of the preceding discussions, the 
global dephasing of the interferometer is obtained by adding 
the contributions of the free propagation segments : 



»Z[AB] 



mm 



%[CD] (22) 



We ignore any dephasing of the probe field associated with 
scattering on the beam splitters. This point will be dis- 
cussed in more details in the next section. 

For simplicity, we consider the interferometer to have 
a rhombic geometry with 2a the aperture angle and £ab 
the length of each side. The latter quantity is related to 
the time of flight tab and the modulus v at of the atomic 
velocity : 



£AB = Wat TAB 



(23) 
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We suppose the interferometer to lie in the plane (x 1 , x 3 ) 
and write the coordinates of its apexes in a manner ex- 
ploiting the symmetry of the rhomb (x 2 = for all points) : 

tD = tab , x D = (^ab cos a, 0, 0) 

t A = —tx> , x A = — x D 

t B = , x B = (0, 0, £ab sin a) 

f c = t B , x c = -x B (24) 

We have used the stationarity to select a specific time. 
Then the reduced velocities are read : 

uab = -^(cosa, 0, sin a) 

c ' 

ubd = — ^-(cosa, 0, — sin a) 



u B d 



c 

UAC 



ucd = uab 



(25) 



We now evaluate the amplitude <fj [AB] corresponding 
to the segment [AB] from equation (J2l|) . We restrict our 
attention to the non relativistic limit v at <cso that : 



2\/2 Huj 



(e 7 .u A b) 2 (■ 



^—iujt-B 



-iujtA. ^ 



(26) 



Using the symmetry of the rhomb, we then deduce the 



amplitude (f k corresponding to the sum over the four seg- 



ments : 



i>Z = 



2i mc 2 

UJ 



uab)' 



5 7 .uac) I sin 



e 1 7 e 3 7 (1 - cos(wtab)) 



a, 



2/; 



(27) 



where X? at is the kinetic energy of the atom measured as a 
frequency. We finally get the apparatus response function : 



At M 



, a 2 , sin 2 (2c*: 



/ 2 (wt A b) 



(28) 



(29) 



where we have introduced the function 



f(x) = 2 (1 - cos(x)) = 4 sin 2 



We now come back to the phase noise variance 
which determines the fringe contrast (|l7]). This expression 
characterizes the decoherence of the atomic interferome- 
ter for an arbitrary gravitational noise spectrum S/Jui], 
for example the spectrum describing the binary confusion 
background |2^,g3|]. As already mentioned, we see that 
this filter may help us to regularize potential infrared di- 
vergence of this spectrum, besides the factor / 2 (wtab) 
which already cuts off low frequencies. 

This discussion can be made more explicit by consider- 
ing the specific case where the spectrum is constant in the 
domain relevant for evaluating the integral ( p"8| ) . Note that 
this is approximately the case for evaluating the effect of 
the binary confusion background for HYPER. As a matter 



of fact, the signal has to be integrated over an averaging 
time ranging between 1 day and 1 month. Such an aver- 
aging time corresponds to the frequency range 1-30/iHz. 
In this frequency range, the binary confusion background 
dominates the other contributions, in particular cosmo- 
logical ones, and it has a quasi-thermal spectrum. Hence, 
the integral (|l8|) can be estimated, at least roughly, by 
replacing the frequency dependent noise spectrum S/,[a;] 
by a constant. Since the binary confusion background de- 
creases at higher frequencies, the result obtained in this 
manner has to be considered as an upper limit for the 
phase noise variance. 

In this simple case, the integral ([l8|) may be deduced 
from the following properties of the function / : 



f 2 (x) = 4f(x)-f(x) 

duj /(ut) _ 1 - e- r \" 



2ir u 2 + r 2 



r 



so that 



2S h Ql t sin 2 (2a 



3 - 4e~ rTAB + e - 2JV AB 

7rT 



(30) 



(31) 



Now, the time of flight tab of the atoms on the segment 
[AB] is of the order of 1.5 s in HYPER |1| and the averag- 
ing time r av is much longer. This entails that the preceding 
formula may be simplified by taking the further limit : 



tab -C r av 



Ttab < 1 



with the result 



-J7 2 t sin 2 (2a)S , ft T A B 



(32) 



(33) 



We note that sin(2a) is the geometrical characteristic of 
the interferometer which determines the difference between 
the two paths and plays the role of the classicality parame- 
ter entering usual expressions of decoherence rates |]30|.[jT| . 

The cutoff r is no longer present in the estimate (B3[) . 
This argument may be laid down in a more general man- 
ner. High sensitivity measurements usually require long 
integration times so that the condition /Vi n t 1 is met 
with ri n t the time of flight of the probe field in the in- 
terferometer. If the noise spectrum StJw] has such a low- 
frequency behaviour that no regularization of the integral 
is needed, then the variance of the phase noise may be 
evaluated by forgotting the filter : 



A<£ 2 ~ [ p-S h [u}A[u] 

J 27T 



(34) 



If, furthermore, the spectrum is nearly flat in the fre- 
quency domain determined by the apparatus function A[u>] , 
then the variance of the phase noise is simply the product 
of the constant value of Sh by the integral of this appara- 
tus function and the latter is obtained from : 



dm f(ujT) 
2ir u 2 



(35) 
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Now we proceed to a numerical evaluation of the de- 
coherence effect. We use the following numbers which cor- 
respond to HYPER with the choice of Cs atoms |51| : 



the transverse direction of the rhomb and it is responsible 
for the beam splitting effect with a change of transverse 
velocity of the atoms : 



m ~ 

"at - 



sin(2cv) = 



133 a.u. ~ 2 x 10 
0.2 m.s" 1 
t 4 x 10 7 Hz 

^trails 



-25 



kg 



0.035 



t'at 



TAB 



*AB 
Wat 



1.5s 



(36) 



ftrans is the transverse velocity communicated to the atoms 
by the beam splitters and mirrors; it will be discussed in 
more details in the next section. Using the noise level in 
the frequency range of interest Sh — 10~ 34 Hz _1 , we finally 
obtain the phase noise variance due to the dephasing of 
the atomic matter-waves in HYPER : 



A<Pi 



10- 



(37) 



Clearly, the decoherence computed in this manner is com- 
pletely negligible, and the fringe contrast is unaffected 
by the direct coupling of gravitational waves to atomic 
matter- waves. 



5 HYPER-like interferometers 

We have already noticed that, in the most sensitive presently 
studied atomic interferometers, optical elements are built 
up on stimulated Raman processes |Q. As a consequence, 
the dephasing seen by the interferometer also picks up the 
gravitation perturbation of the lasers involved in these 
optical elements. In the present section, we give a precise 
evaluation of this 'indirect' effect and show that it largely 
dominates the 'direct' atomic effect studied in the pre- 
ceding section. In order to fix the orders of magnitude, 
we will still consider the numbers corresponding to the 
project HYPER @. 

The use of stimulated Raman processes to build up op- 
tical elements such as beam splitters or mirrors for atomic 
matter waves has been described in a number of papers 
(see for example @||,|0^|,||]). In the stimulated Ra- 
man process, atoms interact with two counter-propagating 
lasers with slightly detuned frequencies fl\ and i?2- The 
detuning is chosen so that the Raman process, absorption 
of one photon in one beam and stimulated emission of 
one photon in the other beam, is resonant with a tran- 
sition between two hyperfine ground states. In contrast, 
the detuning between one-photon interaction and the in- 
termediate excited state is sufficiently large so that spon- 
taneous emission plays a negligible role. This also entails 
that the equivalent duration of the whole Raman process 
is so short that the Raman process can be considered as 
spatially and temporally localized. Besides the transition 
from one ground level to the other, the main effect of the 
Raman process is a momentum transfer between the field 
and atom. This momentum transfer has its direction along 



2ha 



^trans 



phot 



(38) 



where i7 p hot is the nearly common value of the laser fre- 
quencies. 




LASER 

Fig. 2. Schematic description of an atomic interferometer with 
a Mach-Zehnder geometry and beam splitters and mirrors for 
atomic waves built up on stimulated Raman processes : the 
momentum transfer between atoms and photons produces the 
beam splitting; the gravitational perturbation <S> of the dephas- 
ing between the two arms is now picked up not only by matter- 
waves but also by photons. 



The whole interferometer is sketched on figure (§) with 
the atomic and photonic paths now represented. Each mo- 
mentum transfer is accompanied by a change of atomic 
ground state and the atoms are in the same state in the 
output beam as in the input one. This means that the 
amplitudes corresponding to the two arms are able to in- 
terference. However, the gravitational perturbation is now 
registered not only by the atomic matter- waves but also by 
the laser fields. As shown in the following, the 'indirect' 
photonic contribution even dominates the direct atomic 
one. 

We will write the whole dephasing between the two 
arms as : 



(p 



phot 



(39) 



^at is the dephasing picked up by the atomic matter waves 
when they propagate along the linear segments [AB] , [BD] , 
[AC] and [CD] of the interferometer. It has already been 
calculated in the previous section. <P p hot is the gravita- 
tional dephasing of the electromagnetic phases involved 
in the stimulated Raman processes. It can be written as a 
sum over the optical elements : 



phot 



<? P hot [A] - <Z> phot [B] - <2> phot [C] + <£ P hot [D] (40) 



Internal phase factors, corresponding to evolution at the 
different frequencies of the two atomic ground states, do 



8 



Brahim Lamine et al.: Gravitational decoherence of atomic interferometers 



M 



MIRROR 



Fig. 3. Space-time diagram representing the Raman processes 
at the beam splitter B : the photonic lines are the dashed 
lines with a slope unity with reduced space-time variables ct 
and x; the vertical lines represent the motionless macroscopic 
objects constituting the laser sources and the mirrors reflecting 
photons; the nearly vertical heavy line represents the atoms 
moving with a slow velocity v at <C c. The contribution # p hot 
of the photons to the dephasing is due to their perturbation 
by gravitational waves on their paths from the laser sources 
to the atoms. The phases are supposed to be coherent at the 
laser sources, but they are dephased differently on paths L'B 
and LMB. 



not appear in the final expression, due to a proper account 
of energy conservation. But the same energy conservation 
law enforces that the photonic phase ^ p hot is present : the 
change of momentum of the atoms on the beam splitters 
is just equal to the change of momentum of the field state. 

We now evaluate the gravitational dephasing ^ p hot [B] 
of the laser waves involved in the stimulated Raman pro- 
cess at beam splitter B. They are obtained as the effect 
of gravitational waves on the counter-propagating laser 
waves before the latter attain the atoms at point B. We 
use the simplifying assumption that the Raman process 
is instantaneous, so that the dephasing of the two lasers 
has to be evaluated at the same spatio-temporal point. 
The calculation is the same as for standard optical in- 
terferometers Jl5|,[l6| with optical paths described on the 
space-time diagram of figure (||) . As in the preceding sec- 
tion, the macroscopic reference objects, that is the laser 
sources and mirrors reflecting photons, are supposed to be 
at rest initially and to stay at rest, thanks to their large 
mass. 

The coordinates of the point B are fixed as in the pre- 
ceding section. The coordinates of the other apexes are 
displaced with respect to that of B in the direction of the 
axis X3 which is the propagation direction of the lasers 
(see figure ||) : 



tu — tB — T~MB j X M 



CT M B 



*l = *b - (tmb + t lm) , x L = x B - c(t lm - t M b) 
tu = - {t~lm — tmb) , x L ; = Xl 

(41) 

We have introduced the notations tlm for the time of flight 
of photons from the laser source L to the mirror M and 
tmb for the time of flight of photons from the mirror M 



to the atom at B. The reduced velocities have simple ex- 
pressions for photons : 



1 



MB 



-1 



(42) 



If follows from the expressions (0) that all polarization 
factors involved in the evaluation of ^ p hot have the same 

simple form (e r .u) 2 = sin 2 6. 

The dephasing ^ p hot [B] is then deduced from a factor 
(fil [B] representing a decomposition in momentum space : 



phot 



i/>k = P+P- e l "™ B/3 + 



J1"mb/3+ \ _ giw(TLM— tmb)/8- 



P- 



P+ 

ur/B = k^x 1 ^ = ujt B - k.x B 
P± = 1 ± cos 6 



P+P- = sin 2 



(43) 



From now on, we neglect the difference between the two 
laser frequencies fl\ ~ J?2 = ^ P hot- We introduce a nota- 
tion ?7b to represent the phase time of the mode k at the 
point B. We also introduce the notations P± for the factors 
measuring the angle of propagation between gravitational 
and electromagnetic waves. Since these waves propagate 
at the same speed, potential resonances might occur in 
the limit of colinear propagation p+ = or /3_ = 0, as 
shown by the appearance of the denominators in tpk- In 
fact, these resonances do not occur because the numer- 
ators of the fractions and, also, the polarization factor 
P+P- vanish in this limit Q . 

We then perform the same evaluation for the other 
contributions, </^[A], </>^[C] and <fl [D] which are involved 
in the whole photonic dephasing (|40|). Assuming for sim- 
plicity that the geometry is the same for the four beam 
splitters and mirrors, these 4 terms only differ through a 
global phase and their interference leads to : 



A7 _ ^ 



2\/2u 



-lujrjA 



-lujrjB 



-lUJTjC 



-IUJT] D 



(44) 



The phase-times tja , rjc and tjb are defined as t/b from the 
phases of the gravitational mode k at the spacetime points 
corresponding to the passage of the atom at corresponding 
optical elements. 

The magnitude of the photonic dephasing ^ p hot is mainly 
determined by the laser frequency J7 p hot whereas the mag- 
nitude of the atomic dephasing ^ a t was proportional to the 
frequency /2 a t- Since the latter frequency is much smaller 
than the former, it is expected that : 



<2>at < 



phot 



(45) 



We will see at the end of the present calculation that this is 
the case. As a consequence, the correlation between atomic 
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and photonic dephasings will also have a negligible contri- 
bution. It follows that the phase noise variance ( |l8| ) will 
be determined essentially by the photonic contribution : 



(46) 



Note that the two polarizations have the same contribu- 
tion to the present result. 

We use the function ( p9| ) to express the squared am- 
plitudes : 

\&k\ 2 = /(w?7ab) + /Mac) + /Mbd) + /(^Vcd) 

-/Mad) - /Mbc) 
Vab = Vb — Va = tab ( 1 — sin a cos 9 — cos a sin 9 cos ip) 
Vac — Vc — Va = tab (1 + sin a cos 9 — cos a sin cos (p) 

VCD = VAB , VBD = VAC 

Vab =Vb-Va = VAB + VAC 

Vbc =Vc-Vb= Vac - Vab (47) 



and 



\^ k \ 2 = /?+(/?+ - f3-)[f(u((TLM - 7MB)/3-)) 
+/(wtlm/?-) + f(uJT M Bf3+) 
-/(^((TLM - TMB)/3- - TMB/3+)) 
-f(tj{TLM0- + TMB/3+))] + ^/( WT MB(/3+ + /?_)) 

+/3_(/?_-/? + )/(wtmb/?+) (48) 

These equations give the phase noise variance due to pho- 
tons for an arbitrary noise spectrum Sh [oj] ■ 

In order to go further, we perform the same approxi- 
mations as in the preceding section. We consider the case 
of a thermal bath with the noise spectrum Sh constant 
over the frequency domain relevant for the integral (|l8|). 
We focus our attention on the limit of a small r and use 
the integral ( |35|) as well as the further properties of the 
function / : 

f(x)f(y) = 2f(x) + 2f(y) ~ f(x + y) - f(x - y) 

; = Min(|77|, |x|) (49) 



4-7T U} 2 

We deduce from @||) 



A<P 2 phot 4 



^phot^Tphot 

Z 7T 

Tphot = -^(2T(t] AB ) + 2T(r/Ac) - T(r] AD ) - T(r/ B c)) n 

(50) 



Tphot has been defined so that the expression of A<P 2 hot 

has the same form as A<P 2 t in ([33j) ; it is obtained from the 
auxiliary function T which has the same structure as in 



T( V ) =/?+(/?+ -/?_)[Min (\v\,\tlm-t M b\P~)) 
+ Min (IvItlmP-) + Min (\ V \,t M b(3+) 

- Min (|r?|, |(r LM - tmb)/?- - t M b/3+|) 

- Min (|r?|,r L M/3- + t M b/3+)] 
+P\ Min (M,t M b(/3+ + /?-)) 

+/?_(/?_ -/?+) Min (|r,|,TMB/3+)) (51) 

The lengths are of the same order for atomic and pho- 
tonic lines, but the velocity of light c is much larger than 
the atomic velocity u at . Hence the atomic time of flight 
tab is much larger than the photonic ones tmb and tlm- 
It follows that r\ is much larger than the other time pa- 
rameter appearing in the function Min in the preceding 
equation, except for specific gravitational modes propa- 
gating along the segments of the atomic rhomb. Disregard- 
ing these exceptions which have a negligible contribution 
to the spatial mean values, we obtain a good approxima- 
tion for the function T : 

T( V ) ~ /?+(/?+ - /3_)[|t LM - TMB 

- | (r LM - tmb)/?- - tmb/3+1] + + 0-)tmb (52) 

Since the function T no longer depends on the parameter 
rj, the equivalent photonic interaction time r p hot is simply : 

Tphot = j^( T ) n = V t mb 

V = ^</3+(/3+ - 0-) (\x - 1 1/3- - \xf3_ - 2|) 



tlm 
tmb 



(53) 



The numerical factor y = -^si. is a function of the ra- 

a tmb 

tio x = TML which can be obtained through a numerical 

. tmb ° 

integration : 

y(x) 



5ir | 1 _ 3x 2 -3x+l \ f ™ > 1 
2 \2 3z 3 j IOT X ^ 1 



5tt 
12 



for x < 1 



The function ?/(a;) is drawn on figure ||. 

Finally the phase noise variance A<P 2 hot corresponding 
to the photonic lines is essentially determined by the time 
of flight tmb of photons between atoms and mirrors : 



A<P 



phot 



■^phofShJ/TMB 



(55) 



2 7T 

Using the numbers corresponding to HYPER : 

n phot ~ 2 x 10 15 rad.s" 1 , t MB =s 1Q~ 9 s 
tlm^3t M b , Tphot ~2tmb (56) 

we get the following estimation of the phase noise variance 
associated with photonic lines : 



A<P 2 , 

^-**phot ^ 10~ 12 



(57) 
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x 

Fie. 4. Variation of the ratio v = Tphot as a function of the 

° " T M B 

ratio x — : ^ LM - ; y varies from the value tt for x < 1 to the 

TMB ' a 12 — 

value ^ at the limit x 2> 1. 

This indirect contribution due to the gravitational per- 
turbation of electromagnetic waves is much larger than 
the direct effect of gravitational waves on atomic mat- 
ter waves. This result was anticipated from the fact that 
the laser frequency ^7 p hot is much larger than the corre- 
sponding atomic quantity J? a t. It was however necessary 
to perform the whole calculation to reach an unambiguous 
conclusion, because the variance also depends on the times 
of exposition of photons to the gravitational interaction, 
these times tmb and tlm being smaller than tab- 

The final result proves that the decoherence of HYPER- 
like interferometers is dominated by the photonic contri- 
bution. In other words, as far as the coupling to gravi- 
tational waves is concerned, such interferometers essen- 
tially behave as optical detectors with a readout mediated 
by atomic waves. A second important consequence of our 
result is that the photonic contribution to decoherence, 
though much larger than the atomic one, is still com- 
pletely negligible. Hence, the fringe contrast of HYPER 
is unaffected by the scattering of gravitational waves. 

It is worth explaining this result in more details in 
terms of spectra. To this aim, we come back to the ex- 
pression ( |18| ) where the variance is the integral of a phase 
noise spectrum S$ [u>]. We then simplify this expression 
by performing the already discussed approximations : Sh 
is considered as constant and the crossed terms in the 
squared amplitudes 2 are disregarded. Keeping only 
the square terms we get a typical dependence : 

S*[u] ~ S h Q 2 phot l^fi- (58) 

where r is one of the times of flight involved in equation 
(fi"8|). The phase noise spectrum has a magnitude of the 
order of iShi7p tot r and a bandwidth of the order of r _1 
which lead to an integral of the order of Shfip^r. This 
simple estimate fits the result (|55| ) obtained through a 
more rigorous calculation. 

Now the phase noise level may be measured as an 
equivalent vibration noise for the mirrors reflecting the 
lasers. This equivalent noise, written in terms of the posi- 
tion q of a mirror, is approximated as : 

S q [u] ~ S h (cr) 2 ~ KT 34 (m/Vt^y (59) 



We have again used the numbers of HYPER with cr nearly 
equal to lm. This corresponds to a noise level v/Sg ~ 
10 _17 m/v / Hz which is far beyond the vibration noise level 
yJ~S~ q ~ l(r 12 m/VHi which is the target of the HYPER 
instrument. This discussion confirms that the phase noise 
induced by the scattering of gravitational waves is com- 
pletely negligible in HYPER- like interferometers. In par- 
ticular, it shows that this fundamental spacetime noise 
is smaller than the residual phase noise corresponding to 
mechanical vibrations of the mirrors. 



6 Discussion 

It had been suggested that atomic interferometers could 
be sensitive to a decoherence effect stemming from intrin- 
sic spacetime fluctuations. In the present paper, we have 
studied the effect associated with the scattering of gravi- 
tational waves, which we expect to be the dominant source 
of spacetime fluctuations in the frequency domain of inter- 
est for atomic probes. Taking the numbers of the project 
HYPER as an example, we have shown that this effect 
is completely negligible. Essentially, this result has to be 
considered as positive for the project HYPER : as a mat- 
ter of fact, it entails that phase shifts stemming from the 
Lense Thirring effect, the observation of which constitutes 
the main scientific objective of the project will not be 
washed out by the stochastic background of gravitational 
waves. 

The results could be different when considering larger 
sources of spacetime fluctuations or different couplings to 
matter But it is natural to think that the 

gravitational waves, which are predicted to exist in our en- 
vironment by general relativity, are the dominant source 
of spacetime fluctuations in the frequency range of inter- 
est. In any case, the results derived in this paper from well 
established knowledge about gravitational waves and their 
interaction with matter may be used as a reference point 
to which more speculative proposals have to be compared. 

In this concluding section, we discuss a few points 
which may be relevant for a larger class of atomic interfer- 
ometry experiments. In order to discuss the scaling prop- 
erties of the decoherence effect with respect to the main 
relevant parameters, we rewrite the atomic and photonic 
contributions : 

A$l~S h f2 2 at sin 2 (2a) r at 
^phot~^^ hotTphot ( 60 ) 

Sh is the gravitational noise spectrum, supposed to be 
constant over the frequency range of interest, J? a t and 
■f^phot are the kinetic energy of the probes measured as 
a frequency, r at and r p hot represent the times of flight of 
atoms and photons respectively. 

The phase noise variances scale as the times of ex- 
position of the probe to the gravitational perturbation, 
which means that the effect can be understood as result- 
ing from a Brownian-like diffusion due to stochastic fluc- 
tuations of spacetime 0. It is worth noticing that this 
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result is directly linked to the assumption of a fiat noise 
spectrum. A different scaling law for the spectrum, such as 
that predicted by most cosmological models, would nec- 
essarily lead to a different dependence of A& 1 versus the 
time of exposition. 

Now an important point has to be emphasized at this 
stage : contrarily to what could have been expected, the 
key atomic parameter which determines decoherence is not 
the rest energy mc 2 of the atomic probe but rather the ki- 
netic energy mv 2 t /2. For atoms, this makes a significant 
difference in the evaluation of A<P 2 which scales as the 
square of f2 at . For photons, the mass vanishes but the 
kinetic energy J? p hot is merely the frequency. This result 
is directly related to qualitative arguments already pre- 
sented : the phase noise can be expressed through gauge 
invariant expressions and, as a consequence, it may be 
evaluated in a specific gauge, for example the TT gauge; 
in this gauge, gravitational waves have no effect on mas- 
sive particules at rest. Hence, matter-waves corresponding 
to slow atoms are poorly coupled to gravitational waves. 
Clearly, this is not true for the detection of quasi-static 
gravitational field such as the Lense-Thirring effect looked 
for in the HYPER project fl5fl . For such a measurement, 
the sensitivity is effectively determined essentially by the 
rest mass frequency, which explains why atomic interfer- 
ometers may be used as highly sensitive probes of qua- 
sistatic metric effects Jl3|,[53 54 . 



Incidentally, these discussions imply that cold atoms 
are poorly adapted to the detection of gravitational waves. 
Should we aim at observing an effect of stochastic grav- 
itational backgrounds on interferometers, more natural 
strategies would use either optical interferometers or atoms 
with as high a kinetic energy as possible. If the beam split- 
ters and mirrors are built up with stimulated Raman pro- 
cess, this raises the problem that the transverse velocity 
and, therefore, the area of the rhomb decrease when the 
kinetic energy is increased. 

Finally, we want to come back to the qualitative argu- 
ments bearing on Planck units which have been evoked in 
the beginning of this paper. To this aim, we introduce new 
characterizations of the gravitational noise spectrum : 



Sh\u] 



16G i. t r 

-g^g- K B lgw[ 



igg 

~5c5~ 



hujn„ 



(61) 



n gw is the number of gravitons per mode and T gw is an ef- 
fective noise temperature associated with the gravitational 
background. Both quantitites depend on the frequency in 
general. For the binary confusion background considered 
above, we obtain the following values : 



10 

S h ~ 10~ 34 Hz 



Hz < < 10~ 4 Hz 
2" 1 , T gw ~ 10 41 K 



(62) 



The effective noise temperature has an extremely high 
value, even higher than Planck temperature ~ 10 32 K. This 
emphasizes the unconventional character of the noise tem- 
perature T gw from the point of view of thermodynamics. 
In fact, gravitational waves interact so weakly with matter 
that the associated thermalization time is extremely long. 



We also introduce a parameter <9 gw which measures 
the noise temperature as a frequency : 

6) gw =7r ujn gw [Lo} ~ 3 x 10 52 s- 4 (63) 



Using this parameter, we rewrite the phase noise variances 
"301) as : 

A$l t ~ (/2 at i P ) 2 sin 2 (2a) gw r at 



A$ 



at \.""a,t"f ; ^gw 
phot ~ (^phot^p) G'gwTph. 



ot 




(64) 



We still notice the linear dependence of the variances with 
respect to the time of interaction, which is characteristic of 
a Brownian-like diffusion process, as already discussed. We 
also observe the quadratic dependence of the same quan- 
tities in the Planck time tp, which just means that we are 
dealing with effects linear in the Newton constant G. Be- 
sides these two time parameters, the variances depend on 
two frequency parameters fl and <9 gw which measures re- 
spectively the kinetic energy of the probe and the effective 
noise temperature of the gravitational background. 

In the case of the atomic probe, the phase noise vari- 
ance may equivalently be written in terms of the Planck 
mass : 



( mvl t 
V mpc 2 



sin 2 (2a) 6> gw r at 



(65) 



The fraction illustrates the simple scaling argument 
presented in the introduction : for microscopic masses, 
this fraction is much smaller than unity so that the cou- 
pling to gravitational fluctuations tends to become neg- 
ligible. This result has to be contrasted to the fact that 
the scattering of gravitational waves tends to become the 
dominant source of decoherence for macroscopic motions 
[ |37|]3l| . It is however worth acknowledging that the frac- 
tion mt H is only one of the factors which determine the 

phase noise variance A<fr\ t . It is therefore necessary to per- 
form the whole calculation, as we did in the present paper 
for atomic interferometers, before reaching a quantitative 
conclusion about the effect of gravitational fluctuations on 
decoherence. 
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